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We present the first investigation of tlie extrapolation of quenched nucleon magnetic moments in 
quenched chiral effective field theory. We utilize established techniques in finite-range regularisation 
and compare with standard dimensional regularisation methods. Finite-volume corrections to the 
relevant loop integrals are also addressed. Finally, the contributions of dynamical sea quarks to the 
proton moment are estimated using a recently discovered phenomenological link between quenched 
and physical QCD. 



I. INTRODUCTION 

Describing the quark content of nucleon structure in 
terms of QCD is a fundamental aim of modern nu- 
clear physics. As an inherently nonperturbative theory, 
the most rigorous approach to low-energy phenomena in 
QCD is by numerical simulations in lattice gauge the- 
ory. Although restricted at present to the regime of 
quark masses exceeding those realized in nature, recent 
advances in effective field theory (EFT) have made it 
possible to accurately extract the physical nucleon mass 
from QCD m. 

These advances and breakthroughs in chiral effective 
field theory (xEFT) have their origin in the study of a 
range of hadron properties in QCD, including nucleon 
magnetic moments and charge radii Q, H, Hi Hi > the nu- 
cleon sigma commutator [ll B j moments of structure 
functions ^ ^, J^] and the p-meson mass [ll| . 

With the the most detailed studies being on the ex- 
trapolation of the nucleon mass, it has been shown that 
the use of finite-range regularisation (FRR) enables the 
most systematically accurate connection of effective field 
theory and lattice simulation results 0, Il2. Il3| . Mathe- 
matically equivalent to dimensional regularisation (DR) 
to any finite order, FRR chiral effective field theory pro- 
vides a resummation of the chiral expansion with vastly 
improved convergence properties. Central to FRR is the 
presence of a finite energy scale which may be used to 
optimize the convergence properties of the truncated ex- 
pansion. The success of FRR-EFT is highlighted by the 
observation that the higher-order terms of the traditional 
dimensionally regulated expansion, although individually 
large, must sum to zero to describe lattice QCD results. 

Here FRR-xEFT is applied to nucleon magnetic mo- 
ments calculated in lattice QCD. In particular, we inves- 
tigate the modifications required for the extrapolation of 
quenched simulation results. Issues with the formulation 
of EFT on a finite volume are also discussed. We con- 
sider the modifications to chiral loop integrals on a finite 
volume and perform a fixed-volume extrapolation. 

We extend a phenomenological link between quenched 
and dynamical baryon masses [lj| to the case of mag- 
netic moments to estimate the artifacts associated with 



the quenched approximation. We find that the full QCD 
corrections of the quenched magnetic moments are small 
over a wide range of quark mass. 

Finally, the convergence properties of a truncated 
Taylor-series expansion of the FRR results in powers of 
TO^, analogous to that encountered in dimensional regu- 
larisation, are investigated. We illustrate how any moder- 
ate truncation of the series expansion is unable to connect 
with current lattice simulation results. The inclusion of 
higher-order chiral loop corrections associated with the 
A baryon are also considered and found to be small. 



II. MAGNETIC MOMENTS — QUARK MASS 
DEPENDENCE 

In a general construction of an effective field theory for 
low-energy QCD, the expansion of the nucleon's magnetic 
moment about the chiral limit can be written as 
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The term xb-kI-w denotes the leading nonanalytic (LNA) 
chiral correction to the baryon magnetic moment of 
Fig. n Coefficients of low-order nonanalytic contri- 
butions to nucleon properties are determined model- 
independently and are known to high precision phe- 
nomenologically [I^. For example, the LNA coefficient 
to the proton magnetic moment is given by 
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Any extrapolation of lattice QCD simulations must in- 
corporate this knowledge of QCD. 

The analytic terms are unconstrained by chiral symme- 
try, and hence must be determined empirically. Lattice 
QCD provides an ab initio framework to determine these 
parameters from QCD. 

The expression, Eq. |^, is analogous to the chiral ex- 
pansion of the nucleon mass described in Refs. 1, 1^. 
The regulator dependence of the integrals is implicit. 



FIG. 1: Diagram providing the leading non-analytic contribu- 
tion of pions (dashed curve) to the nucleon (solid line) mag- 
netic moment. 



To define the expansion one must define a regularisa- 
tion and renornialisation scheme to remove ultra-violet 
divergences. 

In the text-book approach to elementary field theory 
one commonly uses dimensional regularisation to renor- 
malise the theory. In this case, only the residue of the 
pion pole becomes apparent and the loop contribution is 
therefore given by 
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and the linear divergence of the integral is absorbed into 
an infinite renornialisation of uq . 

Alternatively, ultra-violet divergences can be removed 
by suppressing loop integrals above a finite energy scale. 
The first systematic study of finite-range regularisation 
in xPT '^8'S performed by Donoghue et al. 16]. The 
development of FR R-y PT in the context of lattice QCD 
[1 m in HE in Ii3 has found remarkably improved 
convergence properties of the chiral expansion, providing 
access to a much wider range of range of quark mass than 
the naively regularised theory ^2,] ■ FRR is therefore best 
suited to the problem of chiral extrapolation, where the 
systematic error in the extrapolation of modern lattice 
simulations is of order 1% [l|. 

Using FRR, the loop integral, /^r, in the heavy baryon 
limit can be expressed as 



4 
3^ 



dk 



k'^u^ik) 



(4) 



where uj^ ~ -^//s^ + rn^ and u{k) is the ultra-violet regu- 
lator. 

Expansion of the loop contributions as a power series 
enables one to obtain the renormalised chiral expansion 
parameters, where each of the analytic terms in Eq. (^ 
are renormalised by a finite amount |l2L \\^ . For ex- 
ample, using a dipole regulator, u{k) = (1 -I- fc^/A^)~^, 
Eq. Q) becomes 
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and the Taylor expansion provides 
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Therefore, precisely the same LNA contribution is recov- 
ered 
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with a finite renornialisation of all other terms in the se- 
ries. By varying A, strength in the loop integral may be 
moved to the residual analytic expansion and vise-versa. 
As the moderately large m^ behaviour of the loop in- 
tegral and the residual expansion are radically different, 
varying A provides an opportunity to optimize the con- 
vergence properties of the truncated chiral expansion. 

We show the mathematical equivalence of the renor- 
nialisation prescriptions to a given order. The renor- 
malised expansion in dimensional regularisation is 
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and these renormalised coefficients are recovered from the 
FRR expansion via 
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where the second terms compensate the A dependence of 

In summary, in working to leading order in the chiral 
expansion with a dipole FRR, the quark mass depen- 
dence of nucleon magnetic moments in QCD is 

,B = a^ + aiml + afmt-XB.^^^;^. (10) 

The inclusion of an m^ term here is in anticipation of 
adding NLNA terms in the following. 



III. QUENCHED CONSIDERATIONS 

Here we address the necessary modifications to the chi- 
ral effective field theory for the quenched approximation. 
Vacuum fluctuations of qq-pairs are absent in quenched 
simulations. As a result, the structure of the low-energy 
effective field theory is modified. Meson loop diagrams 
are restricted to only those where the loop is comprised 
of valence quarks having their origin in the interpolating 
fields of the baryon correlation function. This has the 
effect of modifying the effective tt-N coupling constants 
[isL [13 and the corresponding factors xb are changed 
accordingly. 

To summarize the LNA contributions to nucleon mag- 
netic moments in both quenched and dynamical QCD, 
we use the standard notation and define 
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and provide the coefficients, /3^, are in Table HI [TM IT9l | . 
We use U = 93 MeV, D = 0.76 and F = 0.50, {qa = 

D + F). 



TABLE I: Coefficients of the leading pion-loop contributions 
to nucleon magnetic moments in QCD and QQCD. 
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A peculiar feature of the quenched theory is the ap- 
pearance of the flavour-singlet 77'-meson as a light degree 
of freedom. In the absence of vacuum quark loops the rj' 
behaves as a Goldstone boson |23,|21| and must therefore 
be incorporated in the low-energy effective field theory. 

The 77' carries no charge and therefore does not make 
a direct contribution to the magnetic moment of the nu- 
cleon. The leading enhancement of the moment comes 
from Ty'-dressing of the current insertion as in Fig. |21 
Being a flavour singlet, the ry' can propagate through a 
pure gluonic state. This corresponds to a double-hairpin 
in the quark-flow diagram. The vertex correction to the 
magnetic moment induced by this loop produces more 
singular nonanalytic behaviour in the chiral limit than 
the physical theory. 

The double-hairpin diagram has a logarithmic diver- 
gence near the chiral limit. This is a pathological feature 
of the quenched approximation, where the magnetic mo- 
ment tends to infinity near the chiral limit. Physically 
this is not possible as angular momentum quantization 
ensures that the induced magnetic field of the meson- loop 
remains finite, even in the chiral limit ,2Zj . 

The double-hairpin vertex correction diagram has been 
calculated in the graded-symmetry approach to Q^PT by 
Savage in Ref . 18] , and provides the term 
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The superscript (Q) denotes a quenched quantity. The 
coefficient uiq is associated with the double-hairpin ver- 
tex in the rj' propagator. This parameter is related to 
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FIG. 2; Double-hairpin rj' vertex correction to nucleon mag- 
netic moment. 
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FIG. 3: Loop corrections, in units of fiN, evaluated with a 
dipole FRR with A = 0.8 GeV. The dashed and solid curves 
are the pion loop corrections of Fig. in quenched and phys- 
ical QCD respectively. The dash-dot curve is the 77' vertex 
correction of quenched QCD, where the physical magnetic 
moment is used to define the renormalized coupling. 



the physical rj' mass [23, |21| and we choose a value 
rriQ — 0.42 GeV. The loop integral in the heavy-baryon 
limit is given by 



/,,' = — I dk- 



\k) 



oji 



(14) 



The normalization of the integral is such that the LNA 
contribution to this loop is logm^r. 

We also highlight that this contribution is proportional 
to the tree-level moment, /.t^ '°° [l^J. Because of the log- 
arithmic divergence one cannot simply define this to be 
the renormalised moment in the chiral limit. This nec- 
essarily means that one cannot remove the scale depen- 
dence of the coefficient of this chiral log using standard 
methods. We remove the scale dependence by replac- 
ing the tree level coefficient by the renormalised mag- 
netic moment at each quark mass Hg '^'^'^ — /i^ . This 
approximation will be accurate provided Xn' Vs ''^''^77' 
makes only small contributions for 771^ > m^^^^. 

In FigureElwe show the value of the loop contributions 
for varying pion mass. The corrections from the rj' are 
quite small in the region of interest. 

Our expansion in the quenched approximation, analo- 
gous to Eq. (Q), is given by 
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and is used to determine the parameters a^ . The total 
magnetic moment at arbitrary TOtt is then given by 
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FIG. 4: Fit of Eq. I|15^ to quenched lattice data. The solid 
squares (■) illustrate the FLIC fermion results |23 . l23 | and 
the open symbols describe nonperturbatively improved clover 
results (2^], at /3 = 6.0 (V), /3 = 6.2 (A) and (3 = 6.4 (0). 



We note that at this point the decuplet contributions 
have been suppressed as their contributions are higher- 
order in the chiral expansion, when one accounts for the 
octet-decuplet mass splitting realized in the (quenched) 
chiral limit Q . 



IV. EXTRAPOLATION OF LATTICE 
MAGNETIC MOMENTS 

With the expansion of the low-energy EFT determined, 
one needs to fix the values of the unconstrained terms em- 
pirically. Lattice QCD provides information on the quark 
mass dependence of hadronic properties and hence en- 
ables the determination of these free parameters. If the 
data analyzed is within the applicable range of the effec- 
tive field theory then, upon fixing these low-energy con- 
stants, one has an accurate extrapolation to the physical 
regime. 

The electromagnetic form factors of the nucleon have 
recently been studied in simulations of quenched lattice 
QCD |2^ 123, I23. Early simulations of nucleon three- 
point functions have been performed by Leinweber et 
al. in Ref. [2g. Direct lattice calculations of the nu- 
cleon's strangeness form factor have also been carried out 
inRefs. ^^. 

In this study we choose to analyze only the most re- 
cently performed simulations using improved quark ac- 
tions. Results by Gockeler et al. 23] have been obtained 
using the nonperturbatively improved clover fermion ac- 
tion [291 . We also consider recent form factor simula- 
tions by Zanotti et al. 24, 25j u sing the fat-link irrele- 
vant clover (FLIC) quark action [30|- For the purposes 
of this investigation, we select the six most accurate data 
points from the FLIC data set. The precision of the FLIC 
fermion results reflects the use of improved unbiased esti- 
mation techniques j2|| , improved actions and high statis- 



FIG. 5: The dependence of the extrapolated magnetic mo- 
ment determined at the physical pion mass on the dipole regu- 
lator parameter A. Variation in the quenched proton moment 
is similar to the statistical uncertainty of the FLIC simulation 
result at the lightest quark mass considered. 



tics. It has been demonstrated that FRR-^EFT is appli- 
cable up to mj = 0.8 GeV^ [l2| and hence we also choose 
to truncate our data set at this scale of pion mass. 

All previous chiral extrapolations of lattice electromag- 
netic structure have been based on xPT under the as- 
surnption that quenching effects are minimal H, Ij, Hi IE 
[l7l I23L I31I 132 . 1^ . Here we present the first comprehen- 
sive analysis of quenched lattice magnetic moments using 
QxPT for baryon form factors. Brief reports on prelimi- 
nary QxPT extrapolations have recently been presented 
inRefs. 34, 35, 36). 

In Fig. 01 we show fits to quenched lattice data of the 
proton magnetic moment using Eqs. H15(l and (|16|l . In an- 
ticipation of estimating unquenching corrections we have 
adopted the preferred value of A = 0.8 GeV |lj|. The 
logarithmic divergence of Hp in the chiral limit is evident. 

Fig. illustrates the dependence of the extrapolated 
result, evaluated at the physical pion mass, on the choice 
of regulator parameter A. The vertical axis has been fixed 
to that of Fig.Qlin order to display the relevant scale. The 
variation in the quenched proton moment is similar to the 
statistical uncertainty of the FLIC simulation result at 
the lightest quark mass considered. Further reduction of 
the A sensitivity could be obtained by considering higher- 
order terms of the chiral expansion. 



V. EXTRAPOLATION ON A FINITE VOLUME 

The previous section assumed the lattice QCD results 
provide an accurate representation of the results to be 
realized in the infinite-volume limit. However, results 
obtained on a lattice of finite spatial extent will differ 
from those in the infinite volume limit, particularly in 
the chiral regime where the pion Compton wavelength 
can approach the lattice length. Here we extend the for- 



malism of FRR to incorporate finite-volume effects. 

Tlie leading-order finite-volume effect in the chiral ex- 
pansion lies in terms analytic in the small expansion pa- 
f/L, where L is the length of the cubic vol- 
!3j|. Finite- volume corrections also enter through 
the modification of loop integrals. The requirement that 
Greens functions are periodic |37| restricts momentum 
components to the values 
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For the p- wave loop contribution of Fig. ^ where strength 
in the integrand vanishes for fc = as in Eq. Q) , the dom- 
inant effect of discretizing the momenta is to introduce a 
threshold effect jll|, |lj| . Strength in the integrand is not 
sampled until one component of k reaches 27: / L. Since 
chiral physics is dominated by the infrared behaviour of 
loop integrals, the nonanalytic terms of the chiral expan- 
sion exhibit substantial threshold effects. 

To obtain a complete description of the quark mass and 
volume dependence of hadron properties, one must have 
an expansion in both ruq and 1/L |37|- Without data at 
varying lattice volumes it is impossible to determine the 
expansion coefficients of the 1/L contributions. Recent 
calculations have also considered the finite volume cor- 
rections arising from the modification of the leading loop 
integrals to discretised momentum sums pjM l39l |. Im- 
pressive results for a phenomenological description of the 
volume dependence of hadron masses have also recently 
been reported 1^. 

In the absence of lattice QCD results for magnetic mo- 
ments from a variety of lattice volumes, it is not possible 
to rigorously constrain the complete volume-dependent 
expansion. However, it is possible to precisely describe 
the impact of the finite volume on the quark-mass de- 
pendence of magnetic moments on a single, fixed lattice 
volume. 

Evaluating the loop integrals on a finite volume in FRR 
is a rather simple extension. One simply replaces the 
continuum integral by the discrete momentum sum over 
available momenta on a given volume. We formulate the 
finite-volume corrections in the continuum theory and 
therefore continue the momentum sum to infinity. This 
allows the features of finite-range regularization to carry 
over to the finite volume case. For example, Eq. Q be- 
comes 
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The discretized momenta on a cube are given by fc = 
fcminii for fi G Z'^, with the minimum nontrivial momen- 
tum given by fcmin = Stt/L. We note that in obtaining 
the infinite volume form of Eq. Q the angular depen- 
dence of the 3-dimensional integral has been performed 
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FIG. 6: Finite- volume FRR-EFT fit to FLIC fermion results 
for fixed volume. Here the dipole- vertex regulator parameter 
is fixed to A = 0.8 GeV. 



analytically. In general, with an angular dependent inte- 
grand, one must take caution in the naive conversion to 
a spherically symmetric sum as shown in Eq. (|18|l . For 
the integrals used in this paper, we have verified that 
the spherically symmetric and angular dependent sum- 
mations are equivalent. 

Simulations performed with standard Wilson actions 
have large 0(a) errors that can be accounted for in the 
effective field theory. One can introduce new local op- 
erators into the chiral Lagrangian reflecting 0[a) terms 
associated with lattice discretization effects and explicit 
chiral symmetry breaking [4^ 143. li^ . However, there 
has been tremendous success in removing 0(a) errors 
and suppressing ©(a^) errors in lattice simulation results 
through the development of nonperturbatively improved 
actions 30, 44> 451 . These actions display excellent scal- 
ing properties |46l l47l l48l | , providing near continuum re- 
sults at finite lattice spacing. In particular, the FLIC 
fermion simulations, dominating the chiral fits here, are 
performed using an C'(a^)-mean- field improved Luscher- 
Weisz plaquette plus rectangle gauge action |43| and the 
nonperturbatively C'(a)-improved FLIC fermion action 
fSO, 46]. Hence these lattice results already represent an 
excellent approximation to the continuum limit. 

Fig. illustrates the fit of Eq. ((T3|) with the tt and Vj 
loop integrals modified as described in Eq. H18I) for the 
finite volume. The physical volume of the FLIC lattice 
is I^ = (2.56 fm)'^. The chiral properties of the finite- 
volume extrapolation are qualitatively different from the 
infinite volume curve of Fig. 01 In Fig. we show the 
regulator parameter dependence on the finite- volume ex- 
trapolated magnetic moment for A in the range 0.6- 
1.0 GeV. The variation of the extrapolated moment is 
suppressed relative to the infinite volume case, changing 
by 0.1 [ijq . This systematic uncertainty is smaller than 
the statistical uncertainty of the lightest quark mass con- 
sidered here, and could be suppressed further through the 
introduction of higher-order terms in the chiral expansion 
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FIG. 7: The extrapolated magnetic moment on a finite vol- 
ume, V — (2.56 fm)'', for varying regulator parameter A. 



or through the introduction of precise lattice QCD results 
at light quark masses. 



FIG. 8; Correcting the finite- volume quenched approximation 
to the infinite- volume limit of full QCD. The solid curve is the 
finite- volume quenched fit to the data as in Fig. 2] The dotted 
curve provides an estimate of the infinite-volume limit mag- 
netic moment in the quenched approximation. The dashed 
curve shows estimates of the proton magnetic moment in full 
QCD as described in the text. 



VI. ESTIMATING EFFECTS OF DYNAMICAL 
SEA QUARKS 

In a study of baryon masses in quenched and 2-1-1- 
flavour QCD it has been found that the short-distance 
physics of the analytic terms in the residual expansion of 
FRR-EFT are found to be very similar when the chiral- 
loop effects are evaluated with an appropriate FRR |l4l |. 
By identifying the short distance behaviour in QQCD, 
one need only restore the chiral loop effects of QCD to 
obtain an improved estimate of the physical magnetic 
moment. 

In making such an identification it is essential to have a 
consistent method for setting the scale in both quenched 
and dynamical QCD. In particular, one must ensure that 
the procedure is insensitive to chiral physics. The QCD 
Sommer scale |53, based on the static quark potential, 
is insensitive to light quark physics and provides an ideal 
procedure for the scale determination. 

The identification of this phenomenological link be- 
tween quenched and dynamical simulations has been ap- 
plied to FLIC fermion calculations of baryon masses 
|34L |3o|. Upon replacing the chiral loops of QQCD by 
their QCD counterparts the nucleon and Delta are found 
to be in good agreement with experiment. 

By applying the same principle to the calculation of 
magnetic moments in quenched QCD one can obtain im- 
proved estimates of the physical magnetic moment. The 
fit parameters, a^ , are determined by fitting finite- 
volume quenched lattice QCD using Eq. H15|l with dis- 
cretized momenta and a dipole regulator of 0.8 GeV. 
The estimate of the quenching effects are obtained un- 
der the assumption that the bare residual expansion pa- 
rameters are unchanged in infinite-volume QCD when 
A = 0.8 GeV. That is, the fuU QCD result can be de- 



scribed by Eq. (|10(l with the identification a^ = af . 
By fitting with finite-volume FRR-EFT both quenching 
and finite-volume corrections are incorporated in the fi- 
nal estimate. We show the infinite-volume QCD estimate 
of the proton magnetic moment by the dashed curve in 
Fig.El 

In a similar manner, the infinite- volume limit of QQCD 



is estimated by fitting the parameters a. 



S(Q) 



of Eq. ^ 



using finite- volume discretized momenta and a dipole reg- 
ulator of 0.8 GeV in the loop integrals. The correction 
is estimated by Eq. (|16|) calculated with infinite-volume 
continuous momenta in the loop integrals. Figure |H1 illus- 
trates that the finite volume corrections are negligible in 
the regime of the lattice QCD simulation results. 

We emphasize that this result is a phenomenological 
estimate, as the size of the correction is A dependent. 
However an important feature of this approach is that 
the largest finite volume corrections lie in the chiral limit 
as they should. Ultimately, one would like to combine 
the improved convergence properties of FRR-EFT with 
the small l/L expansion such that accurate and model- 
independent determinations of finite volume effects can 
be made. 

The primary feature of Fig. |S1 is that although the 
quenched and physical theory have quite different chi- 
ral structure, the observable effects are rather small. In 
particular, the logarithmic divergence is likely to only be- 
come apparent well below the physical pion mass. Within 
the current formalism of lattice QCD it seems such obser- 
vation would be a formidable task, particularly given the 
large lattice volume required to reveal the 77'contribution. 

The results here, based on the leading chiral correc- 
tions, indicate that proton magnetic moments evaluated 
in quenched simulations give a good approximation to 
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FIG. 9: Truncations of the Taylor expansion of the full QCD 
expression of the magnetic moment, Eq. IjlOfl . to various or- 
ders in rriTr ■ The left most curves display a succession of trun- 
cations from m,r up to m^ . To reach convergence at the 
lightest FLIC data point considered here one requires terms 
up to m^®, while for QCDSF terms to m^ are necessary. The 
dashed curve displays the full QCD expression as displayed 
in Fig. m 



the true theory. The enhancement from the ry'-loop com- 
pensates for the reduction in the standard pion-loop from 
QCD to QQCD. The similarity in the effective curvature 
was also highlighted by Savage [Tg- 



VII. SERIES TRUNCATION AND HIGHER 
ORDER EFFECTS 

It is interesting to explore the convergence properties 
of the series expansion, generated in dimensional regu- 
larisation, truncated to various order in ttitt. Informa- 
tion on the convergence of the truncated series can be 
obtained by studying the Taylor series expansion of the 
FRR dipole regularised form, Eq. ifTUI) . which describes 
lattice results very well. Determination of the conver- 
gence of the series will therefore give insight into the ap- 
plicable range of a dimensionally regularised expansion. 
Similar studies of a truncation of the chiral loop correc- 
tions to the nucleon mass evaluated with a FRR have 
been performed in Refs. |5lll53. l53l| . 

In Fig. El we show the Taylor series expansion of 
Eq. (|10|l truncated at various powers of rn^. This plot 
indicates that in order to accurately reproduce the quark 
mass dependence of the proton magnetic moment up to 
TOtt '-~^ 0.1 GeV then one must keep all terms up to m]^ , 
i.e. m^. If one wishes to reach the quark mass scale 
of modern lattice QCD simulations many more powers 
in m^ are required. In fact, the lightest displayed FLIC 
quark mass requires terms to m^^, similarly the lightest 
QCDSF point needs ml". 

This problem might have been anticipated by the fact 
that at moderately large quark masses the Dirac moment 



FIG. 10: Leading pion-loop contribution to rmcleon magnetic 
moment from the A resonance. 



TABLE II: Coefficients of the leading decuplet, pion-loop con- 
tributions to nucleon magnetic moments in QCD and QQCD 
[13, C = -2D. 

Pba 



Baryon 



Pm 



P 

n 






_ 1^2 

\n2 

6*- 



of the nucleon would be revealed 

eh 



Mp 



2M 



N 



(19) 



Knowing that in this regime, the nucleon mass grows 
linearly with m^ , the moment would require an expansion 
in inverse powers of vn J to describe the data ;2,, 3] . 

We also investigate some of the effects of including 
higher-order terms in the FRR extrapolation of lattice 
data. The leading contribution from the Delta baryon, 
shown in Fig. 1101 is also incorporated into the fit. Equa- 
tion (|15|l becomes 



M^^^ = 



B(Q) 



'(Q), 



S(Q), 



(Q) (Q) J , 



XQ) 



A(Q) 



X^^;/.+x2^^/a^,(2G) 



where, analogous to Eq. Hll() . the corresponding cou- 
plings are given by 






SA ! 



(21) 



and Tabled The loop integral is also modified by the 
fact that the intermediate baryon propagator is non- 
degenerate with the external state. With the mass- 
splitting given by A, Eq. Q) becomes 

^_^r (A + 2^,)fc^ (22) 

Given that the finite-volume corrections are negligible 
in the regime where the lattice QCD results lie, we illus- 
trate the role of the A in FRR-EFT by taking the lat- 
tice results as an accurate representation of the infinite- 
volume limit and evaluate the loop integrals of Eq. H2U|) 
in the infinite volume limit. 



2.5 
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a. 
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FIG. 11: The role of A contributions in the extrapolation 
of the proton magnetic moment. The solid curve displays 
the leading-order QQCD extrapolation of Eq. 11611 . without 
explicit decuplet contributions, as seen in Fig. 2] If Eq. 11511 
is extended to include the pion contributions with decuplet 
baryons in Eq. 12UII . one obtains the dotted curve. 



The fit of Eq. (|20() to lattice results is shown by the dot- 
ted curve of Fig. ^2 The comparison with the leading 
order result (solid curve) shows that the effect of the de- 
cuplet on the extrapolation result is negligible. The new 
non-analytic behaviour introduced by an explicit inclu- 
sion of the A was already approximated well by the ana- 
lytic terms in the expansion at leading order. Although 
the extrapolation shows little sensitivity to the inclusion 
of the A, it will be necessary to explicitly include this 
degree of freedom if one is to extract the low-energy con- 
stants to this order. 

With regard to the estimation of full QCD corrections, 
the inclusion of the decuplet is again a small effect. With 
the A, the QCD estimate of the proton magnetic moment 
at the physical quark mass is increased by 0.05 ^n from 
the leading-order result. Although small for the proton, 
the inclusion of such contributions is found to be impor- 
tant for other baryons of the octet in the extraction of 
the strangeness magnetic moment of the proton |5J, |53 . 

Finally, the prediction of the physical proton magnetic 
moment obtained by including the A contributions and 
compensating for the finite lattice volume is 2.54(30) /iat, 
where the uncertainty is statistical in origin. This result 



agrees well with the experimental value of 2.79 /iat. 



VIII. CONCLUSIONS 

Quenched and physical magnetic moments are in good 
agreement over a large range of pion mass. Although pion 
effects alone are decreased in QQCD, the new behaviour 
introduced by the 77' acts to suppress any difference. 

It is commonly accepted that the dimensionally regu- 
larised expansion of the nucleon mass to 0{p'^) is con- 
vergent to TOtt ~ 300 MeV. We find that to reach sim- 
ilar pion masses in the DR expansion of the magnetic 
moment one would require knowledge of the chiral ex- 
pansion to 0{p^'^). The smooth behaviour of the lattice 
data, together with the series truncations of the FRR ex- 
pansion indicate that although higher-order terms of DR 
can be individually large they effectively sum to zero in 
the region of interest. FRR-EFT provides an effective re- 
summation of the chiral expansion that ensures that the 
slow variation of magnetic moments observed in lattice 
QCD arises naturally in the FRR expansion. 

Finally, by estimating finite volume and quenching ef- 
fects through the leading one-loop contributions of the 
finite-range regularized meson cloud, we obtain an excel- 
lent value for the physical magnetic moment of the pro- 
ton. Combined with the previous success in describing 
the difference between quenched and fuU-QCD nucleon 
and A masses jljj, these results strengthen the argu- 
ment that artifacts of the quenched approximation can 
be accurately corrected using phenomenological meth- 
ods. Ultimately, lattice QCD simulations incorporating 
light dynamical-quark degrees of freedom on a variety of 
lattice volumes are required to achieve the goal of an ah 
initio determination of nucleon properties. 
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